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We perform a geometric study of the equilibrium locus of the flow that 
models the diffusion process over a circular network of cells. We prove 
that when considering the set of all possible values of the parameters, the 
equilibrium locus is a smooth manifold with corners, while for a given value of 
the parameters, it is an embedded smooth and connected curve. For different 
values of the parameters, the curves are all isomorphic. 

Moreover, we show how to build a homotopy between different curves 
obtained for different values of the parameter set. This procedure allows the 
efficient computation of the equilibrium point for each value of some hrst 
integral of the system. This point would have been otherwise difficult to be 
computed for higher dimensions. We illustrate this construction by some 
numerical experiments. 

Eventually, we show that when considering the parameters as inputs, one 
can easily bring the system asymptotically to any equilibrium point in the 
reachable set, which we also easily characterize. 


1 Introduction 


A circular network of cells is a set of cells that are connected in along a ring such that 
the last cell is connected to the first. An instance of this kind of model is the ribosome 
flow model on a ring has been introduced in [0]. 
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Such a model is a parametric dynamical system defined as follows: 


Xi 

X2 




^n—l 

Xn 


XnXn{l - xi) - Aia;i(l - X2) 
AiXi(l - X2) - A2a;2(l - 0:3) 


An— 22 ^n—2(1 Xn—l) An—l^^n—l(l Xn^ 

An—l^^n—1(1 ^'n) An 3 :n(l ^'l) 


The parameters Ai, ■ ■ ■ , An are real strictly positive numbers. They define the degree of 
diffusion between the cells. As shown in [6], if the initial point lies in the hypercube [0,1]”, 
then the system always stays within these limits. Then this system models the occupancy 
levels of a circular chain of n sites (for example a circular DNA), while Ai, ■ • ■ , An are 
transition rates. 

In this work, we focus on the study of the equilibrium locus of this system and in 
particular we study its dependency on the parameters. 

We prove that the set of equilibrium points (Ai, • • ■ , An, ei, • • • , Cn) of the system, 
denoted W, is a smooth sub-manifold with corners of which dimension is n -f 1. 

We also consider the projection onto the first summand tt -.Vn ^ 

(Ai, • • • , An, Cl, ■ ■ ■ , Cn) e-(Ai, • ■ ■ , An) and prove it is a proper surjective submersion. 
Relying on Thom’s First Isotopy theorem j5], we prove that for a given value of the 
parameters, the set of equilibrium points is an embedded smooth curve. Also since the 
base space of the projection is connected, for different values of the parameters, the 
curves are all isomorphic. Considering the fiber over (1, • • ■ ,1) G we show that 

all these curves at connected. 

Moreover, we show how to build a homotopy between different curves obtained for 
different values of the parameter set. This procedure allows the efficient computation 
of the equilibrium point for each value of some first integral of the system. This point 
would have been otherwise difficult to be computed for higher dimensions. We illustrate 
this construction by some numerical experiments. 

This homotopy path tracing method is proven to be valid, as we prove the existence of 
the path the numerical method allows to trace. 

Eventually, we show that when considering the parameters as inputs, one can easily 
bring the system asymptotically to any equilibrium point in the reachable set, which we 
also easily characterize. 


2 Parametric Manifold 


Given the parameters Ai, • • • , A^ G (M+)”, the equilibrium point (ei, • • • , e^) G [0, !]"■ 
must satisfy /(Ai, • • • , A„, ei, • • • , e„) = 0, where: 

^ A„e„(l — Cl) — Aiei(l — 62) ^ 


/(A,e) = 


y Afie^(l 63) ^n) j 
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Here it is worth noting that / has values in and not in since the vanishing 
of the n — 1 hrst components of the vector held implies immediately the vanishing of 
the last one, which is up to the sign merely the sum of the n — 1 hrst components. In 
addition the vanishing of all the n — 1 hrst components of the vector held is easily seen to 
be equivalent to the vanishing of the n — 1 components of the function / dehned above. 

For the sake of simplicity, we have denoted A for Ai, ■ ■ ■ , An and e for ei, ■ ■ ■ , Cn- We 
shall now prove that W = /~^(0) is a smooth manifold with corners. 

The Jacobian matrix of / with respect to A is given by: 


Jx 


-61(1-62) 0 

0 -62(1-63) 

0 


0 6n(l-6i) 

0 6n(l-6i) 

n—1(1 ^n) 6 n(l Cl) 


We want to determine the rank of this matrix. We shall use the following lemma. 
Lemma 1. // 6* = 0 for some i, then 61 = • • • = 6n = 0. If e^ = 1 for some i, then 

61 • • • 6n 1. 


This lemma is algebraically obvious and has a clear physical meaning. If the system is 
saturated in some bin, it is the case as well in all bins. If one bin is empty, so must be 
the other bins. 

Furthermore this lemma has an important mathematical implication. More explicitly 
this says that the only points (A, e) of W for which e lies on the border of the hypercube 
[0,1]” are precisely e = (0, ■ ■ ■ ,0) and e = (1, • ■ ■ ,1). Thus when investigating the 
dependency on e of any claim concerning W, it is will useful to look separately at the 
three following cases: e G (0, l)”,e = (0, ■ • ■ , 0) and e = (1, • ■ ■ ,1). 

For e G (0,1)”, the Jx has full rank, that is n — 1. For e = (0, • • • , 0) or e = (1, • • • , 1), 
Jx = 0. 

On the other hand, the Jacobian matrix of / with respect to e is given by: 


0 A„(l-ei) 

0 A„(l-ei) 

0 A)2 _i (1 Cn) A,2(1 61) -j- Xn—iefi—l 

If e = (0, ■ ■ ■ ,0), then Jg is in row echelon form and no pivot vanishes. Thus in that 
case Jg has full rank. lfe = (I,-- - ,1), a short computation also shows the matrix has 
full rank. 

From this simple argument we conclude that the level set W is indeed a properly 
embedded sub-manifold |3] of (M;^)"" x [0,1]" with corners such that dVn = (M^)"" x 
{(0, • • ■ ,0), (1, • • ■ , 1)} = (M+)’^x{(0, • • • ,0 )}U(M:;)”x{( 1, • • • , l)}. Moreover dim(14) = 

2n — (n — 1) = n -|- 1. 


J. = 


AnCn Ai(l 62) 

XriCn. 


AlCl 0 

-A2(1 —63) A2e2 
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3 The Fibers 


While our primary goal is to study the hber over A for each possible of value of the 
parameters, we start our investigation by considering the hber over a given value of e. 
This will allow in the sequel a deeper understanding of the main subject. 

3.1 Fibers Over e 

Given a point e G [0,1]"", the hber over e is by dehnition the set Ue = {A G : 

/(A,e) = 0}. 

When either e = (0, • • • , 0) or e = (1, • • • ,1), the set Ue is obviously 
Following lemma we consider now a given value e = (ci, ■ ■ ■ , e„) G (0,1)"'. In 
that case the set of Ue is simply of the form with We G since the system 

dehning cxe is linear in A and homogeneous and the rank of Ja is n — 1. In the sequel, we 
shall assume that oj^ is normalized: ||a;e|| = 1, so that it is uniquely dehned. Given the 
geometric form of Ue, we can conclude for two diherent points e and e' both in [0,1]", 
the hbers Ue and either coincide or has empty intersection. In the latter case, their 
closures intersect at the origin. 

3.2 Fibers Over A 

Now let us consider the canonical projection: vr : 14 —)■ (A, e) i—)■ A = (Ai, • • ■ , A„). 

Let us we denote tt the projection from M" x M" onto the first summand M". Then tt is 
the restriction of tt to 14- Also the differential of tt at some point p of 14 is the restriction 
of dn to the tangent space TpiVn) of 14 at p. The projection vr fails to be a local submersion 
at points p where dim(7f(Tp(14))) < n ov equivalently if dim(Tp(14)) nker((i7fp) > 1. This 
would happen if the rank of the Jacobian matrix of / with respect to e at p is strictly 
less than n — 1. Indeed we have the following sequence of facts: 

1 . Tp(W)=ker(d/p), 

2. Tp(14) nker(d7fp) = {(0, • • • - ,m„) | d/p(0, • ■ ■ - ,«„) = 0}, 

3. dim(Tp(14) G ker(d7fp)) > 1 rank( Je(p)) < n — 1. 

Now let us show that rank(Je(p)) < n — 1 never occurs. Recall that: 

0 An(l-ei) 

0 An(l-ei) 

0 Afi_i(l S-n) Afi(l Cl) T Xfi—lCn—l 

Let A the matrix obtained from Je by discarding the hrst column. Then 4 is a square 
(n — 1) X (n — 1) matrix. Let us consider the expansion of its determinant with respect 


Je = 


-AnCn — Ai(l — 62) AiCi 0 

— AnCn —A2(1 —63) A2e2 

XiiCn 0 ■ • • 
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to the first line: 


Aiei[A2e2 • ■ ■ A„_2e„_2(A„(l — ei) + A„_ie„_i)]+ 

(—1)"'A„(1 — ei)[(—A2(1 — 63)) • ■ ■ (—A„_i(l — e„))] = 

Aiei[A2e2 ■ ■ ■ Xn-2Gn-2{^n{^ ~ &i) + A„_ie,i_i)] + 

An(l ~ ei)[A2(l — 63) • • ■ A„_i(l — Cn)] 

When the coordinates of e neither vanish nor are equal to 1 (see lemma [^, it is a sum 
of strictly positive terms. Thus it never vanishes. When e = (0, • ■ ■ , 0) or e = (1, ■ ■ ■ , 1), 
as mentioned in section the matrix Jg has also full rank. This latter fact can also be 
seen directly from the explicit expression of the determinant. 

Eventually, these simple considerations yield the following lemma: 

Lemma 2. The canonical projection tt : W —t is a submersion at every point of 

Vn. 

Note that one could reach the same conclusion by considering the discriminant variety 
of the system |2]. 

Now we turn our attention to the surjectivity of tt and prove the following lemma: 
Lemma 3. The canonical projection tt : Vn —> is surjective. 

Proof. Consider A G Let m < min(Ai, • • • , A„_i). 

Consider r G (0,1). Then we can choose a G (0,1), such that XnCn < mr. Similarly, 
we can choose /3i, • • • , (3n-i ^ (0,1), such that A„a = Aj/3j for alH = 1, • • • , n — 1. 

We shall now prove that we can hnd (ci, • • • , e„) G (0,1)” such that a = e„(l — Ci) 
and f5i = Ciie — Cj+i) for alH = 1, • • • , n — 1, provided we choose r wisely. 

For that purpose, consider the map: cj) : [0,1]” —> [0,1]”, (ei, • • • , e^) h->• (e„(l — 
Cl), ei(l — 62), • • • , e„_i(l — e„)). The map is continuous and dehned on a compact set. 
Consider the subset if = {e G [0,1]” : — !}• Being a closed subset of compact set, 

it is compact as well. Therefore 4>{H) is compact. Now let N or ||.|| be the Euclidean 
norm in M” and let = N o (j. Let p be the minimum of {||x|| : x G 0(Lf)} = ifiH)- 
Of course p G Let u = (ui, • • • ,«„) G [0,1]”, such that V’(u) = p. If p = 0, then 

— Ml) = Mi(l — M2) = ■ ■ ■ = M„,_i(l — Un) = 0. TMs implies that either u = (0, • • • , 0) 
or u = (1, • • ■ ,1). In both cases, n ^ H. Therefore we have p > 0. 

In addition if one denotes by (oi, • • ■ , a„) the standard basis of M”, each segment 
Si = {puj : 0 < p < 1} is mapped by 0 to Sj+i. In order to avoid cumbersome 
notations, we use the following rule: if (fc + l)n > j > kn, Sj (respectively aj) denotes 
in fact Sj-kn (respectively aj-kn)- With this rule granted, let us dehne the triangles: 
/■ = {tpai + (1 - t)pai+i : p,t G [0,1]}. 

Let T C M be the m— dimensional pyramid which basis is H and which top vertex 
is 0. The boundary dT of T is the union dT = /i U ■ ■ ■ U /„ U if. The image of T 
by 0 is the connected set whose boundary is 0(/i) U ■ ■ ■ U 0(/n) U Note that 

0(/*) = {{tp - t{l - t)p^)ai+i + (1 -t)pai+2 : p,t G [0,1]} C fi+i. 

Therefore by the above considerations, 0(T) contains the upper-right quadrant of the 
ball centered at the origin and of radius p, i.e p) fl T. 
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It is now sufficient to choose r such that mr < p. Then (a,/3i, • • • ,/9n-i) ^ -B(0,p) fl 
0(T). Then for the given values of Ai,-- - , A„, there exists e G (0,1)”, such that 
7r(A,e) = A. □ 

Finally we shall prove this last lemma: 

Lemma 4. The canonical projection tt is proper. 

Proof. Let iF be a compact set in We need to prove that 7r~^{K) is compact 

in the subspace topology of W- Since 7t~^{K) = {{X,e) E K x [0,1]” : /(A,e) = 0} = 
{K X [0,1]”) n Vn, it a closed subset of a compact set. Therefore it is compact as well in 
(K X [0,1]”). Thus 7r“^(iF) is also compact in 14- □ 

Stacking together these three results, we get the following theorem: 

Theorem 5. For each A G the fiber of the canonical projection tt is a one¬ 

dimensional smooth manifold. Moreover all the fibers are isomorphic. 

Proof. Since the map vr is a surjective submersion, all the hbers are smooth manifold 
of dimension dim(14) —n = npl — n = l. In addition, since the map is proper, by 
Thom’s First Isotopy Theorem [S], the map is locally a trivial hbration. Since (M^)” is 
connected, the hbers are all isomorphic one each other. □ 

One can wonder if two hbers can coincide or intersect elsewhere than at (0, ■ ■ ■ , 0) and 
(1, ■ ■ ■ , 1). Consider two given points A and A' in If there exists e G (0,1)”, which 

lies on both hbers obver A and A', then relying on the notations of section [3T| we observe 
that A, A' G (Je. Therefore one can conclude with the following lemma. 

Lemma 6. For A and X' in the fibers over X and X' coincide if and only if X and 

A' are collinear, . 

Now we shall prove that all hbers are connected. Given theorem it is sufficient to 
prove that the hber over A = 1 = (I,-- - ,1) is connected. As a matter of fact, in the 
sequel, the vector a\ for a G M will be denoted in bold font as follows a. 

For A G let us simply write tt\ = {e E [0,1]” : (A,e) G 14,7r(A,e) = A} for 

the hber. It is obvious that tti contains the segment [0,1]. We want to prove that tt\ 
is restricted exactly to this segment. For this purpose we shall consider the projection 
0 : 7ri\{0,1} —^ (0, l),e Ci. 

Lemma 7. The map 0 is surjective. 

Proof. Indeed for any value a E (0,1), the point a lies in 7ri\{0, 1} and 0(q;) = a. □ 
Lemma 8. The map is a submersion. 
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Proof. The map 0 fails to be a submersion at points where the kernel of the Jacobian 
matrix Jg (for A = 1) is orthogonal to the vector Oi = [1,0, • • ■ , 0]h 

Since the kernel is a one-dimension real vector space, for this condition to be satished, 
it is necessary that the {n — 1)—last columns of the Jacobian matrix Jg are linearly 
dependent. However, as mentioned above the determinant of this sub-matrix is a sum of 
strictly positive terms. Thus it never vanishes on 7ri\{0,l}. Hence, the map 0 is indeed 
a submersion over 7ri\{0,1}. □ 

Lemma 9. The map 0 is proper. 

Proof. Consider a compact subset K C (0,1). Then K is included in some closed interval 
K C [a, 6] C (0,1). Therefore 0“^(Jr) is a closed subset of 7ri\{0,l}, also included in 
0“^([a,6]). The latter set is closed and bounded in M” and thus is compact. Therefore 
0“^(ii") is a closed subset of a compact and is compact as well in M”. Since tti is closed 
in M", 0“^(it') is also compact in tti. Finally, since 0“^(iF) C 7ri\{0,l}, it is compact in 
the subset topology of the 7ri\{0,l}. □ 

Corollary 10. The fibers of 0 are all isomorphic. 

Proof. Given the threee previous lemmas and by Ehresmann’s theorem [1] , 0 is a locally 
trivial hbration Since (0,1) is connected, all the hbers of 0 are isomorphic. □ 

Therefore, it is sufficient to compute the hber over | to understand the common 
structure of all the hbers. 

So consider that ei = |. Assume that e„ > |. Then given that e„_i = ) ’ have 

Cn-i > Gn- On the one hand, by induction we get a decreasing sequence 62 > 63 > ■ ■ ■ > e^. 
On the other hand, we have 62 = 1 — e„ > e„. Thus Cn < \. Thus we ends up with 
a contradiction. Now let us assume that e„ < |. Then we get an increasing sequence 
02 < 63 < ■ ■ ■ < e„. Together with 62 = 1 — e„, we get this time that > |. All together, 
we get that = Ci = ^. Then it is obvious that Ci = ■ ■ ■ = = |. Hence the hber is 

made of a single point. Therefore all the hbers are singletons. Thus there is a bijective 
submersion from 7ri\{0,1} to (0,1). There the two manifolds are diffeomorphic. Then 
7ri\{0,1} is connected. Stacking this to theorem]^ this yields the following corollary: 

Corollary 11. For each value of X & (M!^)”, the fiber irx is a smooth connected curve. 
Moreover the fiber tti is the segment [0,1]. 

4 Hyperplane Section of A Fiber and Homotopy 
Continuation 


We are now interested in analyzing the intersection between the hber ttx and the hy¬ 
perplane Hs = {e E (0,1)*^ : Obviously, this requires that s G [0,n]. As 

^In this case, we made use of Ehresmann’s theorem, while in theoremj^ we used the Thom first isotopy 
theorem, which is a generalization of it. In the latter case the use of the more general setting was 
necessary since F„ is manifold with corners, while here 7ri\{0,1} is a mere manifold and the classical 
setting of Ehresmann is enough. 
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mentioned in j6], for each s and each A, there is a single equilibrium point in Hg. Moreover 
starting from any point in Hg, the system converges toward this equilibrium point. 

Here we shall prove that indeed the intersection Vg^n = Hg (iVn is a smooth manifold 
isomorphic to (M”)"'. More precisely, we have the following proposition. 

Proposition 12. For s G (0,n), consider the intersection Vg^n = Hg (1 Vn- Then Vg^n 
is a n—dimensional smooth variety and the restriction of the projection vr to Vg^n is a 
diffeomorphism. 


Proof. First we shall prove that Vg^n is indeed a smooth manifold. The set Vg^n is defined 
as /“^(O) n hj^{0) for hg{\,e) = ~ ^ ^ (0,■«-) and any n > 1, Vg^n 7^ 0- 

Indeed for a given value of A, the fiber a\ is a connected curve that joins the two extreme 
points (0, ■ ■ ■ ,0) and (1, • • • ,1), which are located each in each half space that Hg fl M"' 
defines. 

Consider the function g{\,e) = (/(A,e), hg(\,e)) G MP and its Jacobian matrix at a 
point (A,e) G Vg^n- 


W = 


— AnCn — Ai(l — 62) Aiei 0 

— AnCn —A2(1 — 63) A2e2 


0 A„(l-ei) 

0 An(l - ei) 


AnCn 

1 


0 

1 


0 A^_i(l Cfi) An,(l Cl) -j- Afi—ien—i 

1 1 1 


Above in section |3i2l we have introduced the matrix A which is the matrix obtained by 
discarding the first column and the last row of W. We have shown that A is invertible for 
any value of (A,e) G 14- Let v be the first column of W without the last element. Since 
rank(y4) = n — 1, v is a linear combination of the columns of A. Moreover v is collinear 
with no column of A. Therefore det(A — V) ^ 0 where V is the (n — 1) x (n — 1) matrix 
which columns are all identically v. As easily seen det(lF) = det(A — V) 4 0. Thus W is 
invertible for all (A,e) G IA,n- Therefore IA,n is indeed a smooth variety of dimension n. 

Now consider the restriction of the projection vr : I4 —)■ to Vg^n- It is a proper 

as 71 is. It is surjective, since as argued just above each fiber of vr meets Hg and vr is 
surjective. It is a submersion because tt is a submersion and the kernel of dir at (A,e) is 
precisely the tangent space to the fiber and Vg^n is transverse to each fiber in I4. Indeed 
the tangent space to Vg^n is precisely the direction of Hg, which is transverse to I4 in 
since rank(IT) = n. 

Eventually, we can conclude that 7i\Vs n is n proper surjective submersion. Again by 
Erhesmann’s theorem, its fibers are all isomorphic since is connected. But over 

(1, ■ ■ ■ ,1) the fiber is a singleton, since the fiber of vr is a segment that intersects the 
hyperplane Hg at a single point. As a consequence is a bijective submersion and is 
therefore a diffeomorphism. □ 

We shall rely on this proposition to prove the existence of the construction we introduce 
now. 

In the simple case of tti, we have immediately fl tti = {-}. Here (as in the sequel) 

= 1 - 
Vn’ 


we have used the rule introduced in section 3.2 


so that - 

ti 


-j G 


In this 






section, we shall show how to hnd the equilibrium point tt^ fl Hg for any admissible value 
of the parameters A. 

Consider a straight segment line in (M^)” joining Ai = 1 to another point A 2 ; X(t) = 
(1 — t)Xi + 1 X 2 for 0 < t < t. This yields an homotopy between the hber n Hg and 
the fiber ttaj H Hg\ 

H{e, t) = (1 - t)g{Xi,e) + tg{X 2 , e), 
where g{X,e) = [/(A,e), e* - s]* e 

As we have seen solving g{Xi,e) = 0 is extremely easy. For a given s G (0,n), starting 
from the solution {-}, we can follow a path embedded in Hg that ends at the solution of 
g{X 2 ,e) = 0. Consider 'j{t) = [xi(t), ■ ■ ■ ,Xn{t)Y such a path. By construction, we want 
it to satisfy for all t G [0,1]: H{'j{t),t) = 0. This yields the following differential system: 


n 


E 


m 

dxj 




dxj 

dtt 


(t) + 


dH, 

dt 


( 7 (f), t) = 0 , 


where Hi,-- - ,Hn+i are the components of H. This is a special case of homotopy 
continuation methods jl]. This system yields a Newton like method for tracing the path 
7 (f), since we have: 




where is the Jacobian matrix of H computed for the n first coordinate at the 

point 7 (t). One can easily write this differential equation explicitly. 

The existence of a solution to this non-linear equation on [0,1] is granted as follows. 
The function g is linear with respect to A. Thus if( 7 (f),f) = g{{l—t)Xi+tX 2 ,'y{t)) = 0. 
Therefore 7 (f) is always an equilibrium point for each value of t. Hence the curve 
f H->• ((1 — t)Ai -|- 1 X 2 , ■yit)) is a lift of the segment joining Ai to A 2 on the variety Vg^n- 


Such a lift exists and is unique by proposition 12 


In our case, the system is particularly simple since ii is a quadratic function of y{t). 
The procedure can be easily implemented. We have implemented the path tracing 
procedure in Python. In the hgure 4T, we show the homotopy path generated for a 
system with 3 bins. 

We have also run the homotopy path tracing procedure with n = 100. The procedure 
allows to hnd efficiently the right equilibrium point in the second hber, which would have 
been otherwise difficult to determine. 


5 Control 


As shown in section 3.1, for a given value of a desirable equilibrium e = (ei, • • • , e„), 
one can hnd easily the values parameter set A that would lead the system to this point, 
provided the initial state a;( 0 ) = (a;i( 0 ), • • • ,x„( 0 )) satishes = Y'i ^i- 

Hence the reachable set for a given initial state vector x(0) is obviously the intersection 
between the cube [ 0 , 1 ]” and the hyperplane dehned by Yi open-loop 

command that would asymptotically lead to an equilibrium point e lying in the reachable 
set can be dehned by taking any value of the parameter set from a^. 
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I — Homotopy Path] 



Figure 4.1: The homotopy path generated with the following parameters Ai = (1,1,1), A 2 = 
(1.39328599,8.30098374,3.98355604) and s = 1. The first equilibrium point Eq = 
(1/3,1/3,1/3) is represented by a green ball, while the final equilibrium point 
El = (0.53112814,0.1203633,0.34850856) is rendered as a red triangle. Of course 
all over the path the constraint ei + 62 + 63 = s holds. 


6 Conclusion 


We have studied the equilibrium locus of the flow generated by a circular network of 
an arbitrary number of cells. Such a model appears for example in the ribosome flow 
model on a ring. We have proven that for every value of the parameter vector A, the 
equilibrium locus is a smooth curve. Moreover, we have shown how to effectively compute 
the equilibrium point for a given value of A and a given value of hrst integral made by 
the sum of system coordinates. 

We have proven that when considering the set of all possible values of the parameters, 
the equilibrium locus is a smooth manifold with boundaries, while for a given value of 
the parameters, it is an embedded smooth and connected curve. For different values of 
the parameters, the curves are all isomorphic. 

Moreover, we have shown how to build a homotopy between different curves obtained 
for different values of the parameter set. This procedure allows the efficient computation 
of the equilibrium point for each value of some hrst integral of the system. This point 
would have been otherwise difficult to be computed for higher dimensions. We have also 
performed some numerical experiments that illustrate this construction. 
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Finally, we considered the control of such a system in open-loop. We showed how the 
system can be driven asymptotically to any reachable equilibrium, when the parameter 
set can be viewed as a vectorial input. 
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